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ABSTRACT 



Aims. We measure the thermodynamic properties of cluster outer regions to provide constraints on the processes that 
rule the formation of large scale structures. 

Methods. We derived the thermodynamic properties of the intracluster gas (temperature, entropy) by combining the 
Sunyaev-Zel'dovich thermal pressure from Planck and the X-ray gas density from ROSAT. This method allowed us to 
reconstruct for the first time temperature and entropy profiles out to the virial radius and beyond in a large sample of 
objects. 

Results. At variance with several recent Suzaku studies, we find that the entropy rises steadily with radius, albeit at at a 
somewhat lower rate than predicted by self-similar expectations. We note significant differences between relaxed, cool- 
core systems and unrelaxed clusters in the outer regions. Relaxed systems appear to follow the self-similar expectations 
more closely than perturbed objects. 

Conclusions. Our results indicate that the well-known entropy excess observed in cluster cores extends well beyond the 
central regions. When correcting for the gas depletion, the observed entropy profiles agree with the prediction from 
gravitational collapse only, especially for cool-core clusters. 
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1. Introduction 

The outskirts of galaxy clusters are characterized by ongo- 
ing accretion processes and significant departures from viri- 
alization and hydrostatic equilibrium. Still, at the present 
epoch they host the transition between regular large- 
scale accretions and the injection of shocks and turbu- 
lent motions, and they are the sites where the tran- 
sition between large-scale infall motions and the dense 
and hot ICM occurs. Departures from hydrost atic equilib - 
rium i n the form of turbu lence motion (e.g., iDolag et~al.l 
120051: IVazza et al.l l2011a^ and non-thermal sources of 
energy such as cosm ic rays (e.g., IVazza et al.l 2012a : 
Pfrommer e t al.l 20071) an d magnetic fields (jPolag et al.l 
1999; Briigg eiTet al.l l2005h are expected to play a more 
important role in these regions than in cluster cores, and 
the infall of smaller structures along large-scale filaments 
may cause the material in these regions to be clumpy 



is required. However, because of the very low surface 
brightness of the X-ray signal, providing reliable tem- 
erature measurenients i n these regions is challenging 
Ettori fc Molendil l201lt lEckert et all l201lt) . Recently, 
the low background of the Suzaku satellite allowed for 
the measurement of temperature profiles out to i?2n[fl 
handful of systems (e.g., iReiprich et al.l 12009 : 



I 



2010|;_ 



Simionescu et al.l 
20121) . 



m a nandtul q i sy stems I e.g.. IKeipricn et al.l izuua 
Bautz et al. 20091: iKaw ahar ada et alT 20101: iHoshino et al.l 



12011 lAkamatsu et al.l 12011 



(jNagai fc Lau[ 20111: Simiones cu et al 



ric (Vazza e t al.l l2oTl b: Eck crt et al 



2m) and asymmet- 
1)12). These effects 



may bias the meas urements of cluster masses using X- 
ray and SZ proxies (iRasia et all 12004 iNagai et al.l[2007bt 



iPiffaretti fc Valdarninill2008HBurns et al. 11201 Ot) , thus limit- 
ing the use of gal axy clusters as high-precision cosmological 
probes (I Allen et al. 2011). 

For all these reasons, a deeper understanding of 
the state of the intracluster gas in cluster outskirts 



Humphrey et al.l 120121 ). However, these measurements 
are often probing a small fraction of the total extent 
of these objects around the virial ra dius, and they ma y 
be subject to systematic problems (jEckert et al.l 1201 ll ). 
Hence, the general behavior of the temperature profiles 
beyond i?5oo is still largely unknown. Conversely, while 
direct temperature measurements require high-quality 
X-ray spectra and an excellent knowledge of the systematic 
uncertainties, estimates of the (projected) gas density 
can be easily obtained from the flux in the soft band 
(0.5-2 keV), provided that the temperature does not fall 
below ^ 1.5 keV. In this respect, the large field-of-view 
of i?OS'^r/PSPC has allowed for the reconstruction of 
gas density pro files out to the vir i al radius in relatively 
large samples (jEckert et all 120121: IVikhhnin eraD 119991: 



^ For a given overdensity A, 
MA/(4/37ri?i) = Apc 



Ra is the radius for which 
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iNeiimannl [2001 . providing a consistent indication of a 
steepening of the gas density profiles beyond i?5oo- 

While the X-ray signal is proportional to the square 
of the gas density, and thus decreases sharply in 
clust er outskirts, the SZ effect (jSunvaev fc ZeldovichI 
|1972| ) depends on the integrated pressure along the 
line of sight, which decreases more gently vifith radius 
(jArnaud et al.l [2blCl[ ). New-generation SZ instruments are 
now routinely detecti ng galaxy clusters, sometimes out 
to la r pe radii (e.g., Plagge et al. l2010t Morandi et al.l 



[20TI lAMI Consortium et all 12011 ISavers et al.l I2012D 

Recently, measurements of the SZ effect in a sample of 
62 galaxy clusters with the Planc k satellite were reported 
(jPlanck Collaboration et al] l2012l hereafter P12), detect- 
ing a stacked signal as far as ~ 3i?5oo from the cluster cen- 
ter. Planck therefore provides high-quality pressure profiles 
in cluster outskirts, which can be combined with X-ray in- 
formation to get some important insight on the state of the 
intracluster medium (ICM) beyond i?5oo- 

In this paper, we demonstrate that by combining the 
average SZ pressure profiles from P12 wi th the aver- 
age R OSAT gas density profiles presented in lEckert et al.l 
(|2012l hereafter E12), it is possible to accurately recover 
the temperature and entropy profiles in galaxy clusters, 
bypassing the use of X-ray spectroscopic information (see 
also Amcglio et al. 2007, 2009). This allows us to probe the 
importance of additional effects (non-thermal pressure sup- 
port, gas clumping, asymmetry) on the quantities recovered 
by assuming that the gas is in hydrostatic equilibrium. In a 
companion paper (hereafter Paper II), we use this informa- 
tion to infer the average hydrostatic mass and gas fraction 
profiles. Similar results are also presented for a subset of 18 
individual clusters that are in common between the samples 
of P12 and E12. 

Throughout the paper, we assume a ACDM cosmology 
with n„i = 0.3, f7fc = 0.05, and Hq = 70 km s^^ Mpc'K 

2. Method 

2.1. Basic formalism 

For this work, the two main observables are the electron 
pressure Pe{r) from SZ measurements and the number den- 
sity of electrons ne(r) from X-ray observations. The gas 
temperature can be obtained from these two quantities as 



kT{r) = 



Heir) ' 



(1) 



with He ^ 1.21n^f. Similarly, the entropy K can be written 
as 



K{r) = kT{r)n,{r)-^^^ = P,ir)n,{ry'^\ 



(2) 



From the self-simil ar model (e.g., iBrvan fc Normanlll998t 
lArnaud et al.ll200"2l ). these quantities are expected to scale 
like 



T - M2/3/i(2)2/3 

K ~ m2/3;i(z)-2/3^ 



(3) 
(4) 



where h(z) = ■\/rim(l -I- z)^ + ^Ia, until the accretion ra- 
dius (~ 1.5 — 2 times the cluster's virial radius), and thus 
we rescale the profiles by these factors, such that in a purely 
self-similar case all profiles should match. 



The radial dependence of the pressure is assumed to 
follow a generalized Navarro-Frcnk- White (GNFW) profile 
(Navarr o ct al...l997; .Nagai et al.„2007ai) . 



Po 



Pjr) _ 

^^500 (c50oa;)^(l + (c5ooa;)")('^-'^'/" '' 



(5) 



where x ~ ?'/^500i and P500 is the expected self-similar 
pressure at i?5oo (see Appendix A of lArnaud et al.l [20101 
for details). 



P^ 



500 



1.65xl0"3/i(z)^/3 



500 



3 X 10i4Af, 



o 



2/3 



/i?o keV cm"^(6) 



In addition, the pressure profiles are rescaled by the quan- 
tity 



500 



3 X 1014M, 



0.12 



(7) 



which acco unts for the de pendence of the g as fraction on 



mass (e.g.. lArnaud fc E~ ard 19991; lEttori fc Fabian 1991 



IVikhlinin et al.ll2006D . In Fig. [CT] and [C2] of the online 



version we show the best-fit GNFW profiles compared to 
the Planck data. The best-fit parameters are found to be 
always consistent with the results given in P12. 

To describe the gas density profiles, we use a simpli- 
fied ve rsion of the functional form given by IVikhlinin et al.l 

(Hoi, 



rigasir) = no 



(l + r2/r2)3/3/2 {l + r-f/riyh 
R 



(l-Hr2/r22)3/3/2 J ' 



(8) 



where Ugas = fie -\-n-H- This function is projected along the 
line-of-sight £ assuming spherical symmetry, and is fit to 
the observed ROSAT emission-measure (EM) profiles. 



EM{r) = 1.21 / n^(w) de, with uj"^ ^ + f 



(9) 



The para meters in Eg. [51 are str ongly degenerate, therefore 
following IVikhlinin et al.l (|2006l ). we fix the parameter 7 to 
the value of 3. We stress that this functional form can allow 
for both a steepening and a flattening of the density pro- 
files, through the sign of the e parameter. This function is 
found to adequately fit the observed emission-measure pro- 
files in aU cases {Xred — !)■ -f'or comparison, we also used the 
density pro files extrac ted th rough direct depro jection as in 
E12 (Kr iss et al.lll983D . The lKriss et~all (|1983D geometrical 
deprojection technique assumes that the outermost density 
bin is unaffected by projection effects, which leads to fiatter 
profiles at the edge. To alleviate this problem, we applied a 
correc tion for the edge effect as explained by McLaughlin 
(|l999f ): iBuotel (|2000f ). After applying this correction^ the 
parametric and deprojected density profiles are found to 
match very well, showing differences that usually do not 
exceed 10% at all radii. The comparison between the de- 
projected and parametric density profiles can be found for 
each individual cluster in Fig. IC.3I and IC.4I of the online 
material. 
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Fig. 1. Thermal pressure (left, using the best-fit parameters from P12) and gas number density profiles (right, from E12) 
for the 18 objects in common between the two samples. The profiles are drawn using the functional forms defined in Eqs. 
[5]and|8l respectively. In both cases, blue lines indicate CC systems, while NCC profiles are shown in red. 



In Fig. [T] we show the pressure (left) and density pro- 
files (right) for a subset of 18 clusters which are in com- 
mon between P12 and E12. Based on their central entropy 
(|Cavagnolo et al.ll2009D . six of these clusters are classified 
as cool-core (CC, Kq < 30 keV cm^), while the remain- 
ing 12 are non-cool core (NCC). The list of these clus- 
ters is given in Table lA.ll We note that in our sample, 
a detection of the emission beyond i?200 a-t a significance 
of at least 2a is achieved for only three CC clusters, so 
the average CC profiles may not be representative in this 
radial range. To rescal e the density profiles, we u s ed th e 
values of -R500 given in iPlanck Collaboration et all ()201lD . 
For the general pro perties of th ese clusters, we refer to 
iPlanck Collaboration et al.l (|201l[) and E12. 

2.2. Monte-Carlo Markov chain estimation of errors and 
confidence intervals 

Providing accurate and meaningful error estimates and con- 
fidence intervals is crucial to probe the state of the ICM 
in cluster outskirts, where both the X-ray and SZ signal 
are weak. We adopt here a Monte-Carlo Markov chain 
(MCMC) technique to compute the confidence envelope of 
the functional forms describing the gas density (Eg. [gt and 
the pressure (Eq.[5]). MCMC (e.g.. lGilks etaLlll996[ ) is the 
method of choice to probe the error distribution of parame- 
ters in fitting multi-dimensional functions. It is particularly 
suitable to cases where the parameter space is complex and 
parameters are strongly correlated. MCMC also allows one 
to make any statistical inference on the basis of the distri- 
butions it has sampled. 

We developed an implementation of MCMC based on 
the Metropolis-Hastings algorithm (formally Metropolis al- 
gorithm, because the Hastings ratios were always 1). Chains 
were calculated for all parameters of the pressure and den- 
sity profiles for all 18 objects and for the average profiles. 
Gaussian proposals were used, taking into account the one- 
dimensional errors on each parameter, which were rescaled 
to reach a rejection ratio between 0.5 and 0.8 for each chain. 
Each of the 38 chains contains 10^ steps. No burn-in phase 
was used, because the chains started from the best fit; how- 
ever, we checked that discarding the first 10^ steps had no 



effect on the chain distributions. Since we are mostly in- 
terested in the cluster external regions, the procedure was 
applied only on the data beyond 0.2i?5oo. The parameters 
affecting only the central regions (< 0.1i?5oo) were fixed 
to their best-fit values during the MCMC procedure; these 
parameters are 7 for the pressure profiles, which was fixed 
to the value of 0.31 as in P12, and and /3 for the density 
profile. Since for both functional forms the parameters are 
strongly correlated, this allows one to explore the parame- 
ters relevant to the outer regions in more detail. In the case 
of the average NCC pressure profile, fixing the inner slope 7 
to the value of 0.31 leads to an underestimation of the data 
points by 10-15% at all radii beyond 0.2i?5oo. Since our fo- 
cus for this work is on the outer regions, in this specific case 
we allowed this parameter to vary during the chain, which 
gives a better description of the data beyond 0.2i?5oo. We 
checked the convergence of the chains by running a few of 
them for 2 x 10^ steps and comparing the distributions of 
the short and long chains. 

Uncertainties in the fit parameters can be straightfor- 
wardly propagated in the determination of the temperature 
and entropy profiles using the chains, allowing one to de- 
termine the confidence intervals on these profiles. For the 
average profiles and for each object, the temperature and 
entropy profiles were constructed for each step of the corre- 
sponding chains. The confidence interval at a given radius 
is then given by the distribution of the profile values at this 
radius. For example, Fig.[5]shows the error distributions for 
the temperature (left) and entropy (right) at i?200 from the 
average pressure and EM profiles built using MCMC. 



3. Results 

3.1. Validation of the method 

To validate our method for the reconstruction of gas tem- 
peratur es, we used temperature measurements from the lit- 
erature (iLeccardi fc Molendi|[2008bt ICavagnolo et al.ll2009l: 



'Sn owden et al.l 2008^ and compared our temperature pro- 
files reconstructed through Eq. [1] with the measured tem- 
peratures. Given that most of these measurements are lim- 
ited to the central regions (< 0.5i?5oo), a similar check 
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Fig. 2. Posterior distributions for the reconstructed tem- 
perature (left) and entropy (right) at i?200 from the av- 
erage pressure (from P12) and density profile (from E12) 
after 10^ MCMC steps. 



cannot be performed for the outer regions. In Appendix 
[Bl we show the temperatures for all 18 objects compared 
with the available data. Outside the innermost ~ 1 ar- 
cmin, good agreement is found between the measured and 
reconstructed temperatures. Our method is typically able 
to reproduce the observed temperatures with an accuracy 
of < 10%. In the innermost regions {R < 0.1i?5oo), angular 
resolution effects may play a role, and thus we excluded 
these regions from our analysis. 

We also tested whether the average profiles extracted 
from the P12 and E12 samples agree with the individual 
objects. Indeed, when average quantities are extracted from 
different samples, selection effects may play a role, in par- 
ticular in the central regions where large cluster-to-cluster 
variations are found. In Fig.[3]we show the temperature pro- 
file obtained by computing the median of the temperature 
profiles from each individual object and compare it with 
the temperature profiles from the universal pressure profile 
(P12) and the median gas density profile (E12). The uncer- 
tainty in the median was computed at all radii through 1000 
bootstrap reshuffling of the populations. At each radius, we 
randomly reselected a set of 18 data points with repetition, 
and computed the median of the shuffled dataset. The me- 
dian value and its error are then given by the median and 
the root mean square deviation of the distribution obtained 
in this way. We note that this error estimation method is 
sensitive to the scatter of the data points around the mean 
value, i.e., a larger scatter (statistical or intrinsic) leads to 
a larger error on the median value. 

In both cases, we show the temperature profiles ob- 
tained from the best-fit parametric function to the EM 
profile (Eq. [5]) with its associated MCMC error envelope 
(Icr confidence level), and by interpolating the deprojected 
density profile using a cubic spline method. The difference 
between these two curves thus highlights the uncertainties 
associated with the different deprojection techniques. All 
methods are found to give consistent results around i?200, 
where the scatter of the values obtained with the different 
methods is consistent with zero. In the outer regions, the 
uncertainties on the median associated with the parametric 
forms are smaller than with the deprojection. This effect is 
caused by the larger scatter of the data points in the latter. 
Indeed, in the parametric approach the shape of the profiles 
is forced to be regular, which reduces the cluster-to-cluster 
scatter and thus the errors on the median (see above). For 
this reason, the profiles obtained through deprojection stop 
at lower radii. 

In the central regions (< 0.3i?5oo), the temperature ex- 
tracted from the average samples falls below the one ob- 



'200 




Fig. 3. Average self-similar scaled median temperature pro- 
file obtained over 18 individual objects compared to the 
profile derived by combining the universal pressure profile 
from P12 and the median gas density profile from El 2. In 
both cases, we show the comparison between the profile ob- 
tained using the functional form for the gas density (Eq.[8l), 
with MCMC error bars, and using non-parametric depro- 
jection. 

tained for the individual systems. This probably indicates 
the presence of different selection effects in the two samples. 
Namely, the X-ray sample is likely constituted of objects 
with a more peaked density profile in the central regions 
compared to the SZ sample, which leads to a decline in the 
temperature. In the outer regions, the temperature profiles 
reconstructed from the average profiles appear to be some- 
what steeper than the average of the 18 individual systems. 
This may once again reflect a selection bias. In all cases, the 
profiles extracted from the 18 systems should be preferred, 
since their selection is much better controlled. 

For the remainder of the paper, we present only the 
results obtained with the parametric form for better read- 
ability. However, as shown in Fig. |21 the profiles obtained 
through geometrical deprojection are found to agree very 
well with those obtained with the parametric form, since 
at all radii the differences are smaller than the statistical 
uncertainties. This conclusion is also valid for the other de- 
rived quantities. 

Another possible source of uncertainty in our 
temperature determination is the conversion between 
i?05'^T/PSPC count rate and EM in the low-temperature 
regime. Indeed, if the temperature in the outer regions falls 
below ^1.5 keV, this factor depends strongly on the tem- 
perature, while at higher temperatures the dependence is 
mild. To take this effect into account, we used the Planck 
pressure profiles and Eq. [T] to obtain an estimate of the 
temperature at each radius assuming temperature profiles 
from the literature, as in E12. Then we used this tempera- 
ture to recompute the conversion factors at each radii and 
extract again the density profiles and the derived quanti- 
ties. This procedure leads to small differences (< 10%) in 
the conversion factors, which converts into differences of less 
than 3% in the reconstructed density. On the other hand, 
line cooling becomes important especially at low X-ray 
temperatures, and thus the assumed metallicity can have 
an impact on the conversion factors. Several works (e . g . , 
iDe Grandi fc Molendil[200ll: iLeccardi fc Molendill2008aD in- 
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Fig. 4. Average self-similar scaled temperature profiles for 
the CC (blue) and NCC (red) populations separately. The 
shaded areas show the results obtained by combining the 
average CC/NCC pressure profiles from P12 with the cor- 
responding profile from E12, while the data points show the 
median of the individual systems in common between the 
two samples. 



dicate that outside the central regions the metallicity pro- 
files are flat with a mean abundance of ~ 0.25Zq; however, 
in these works the metallicity profiles do not extend be- 
yond ^ 0.5i?5oo. Because of the poor spectral resolution of 
ROSAT /PSPC, estimating the behavior of the metallicity 
profiles in the outer regions is beyond the scope of this pa- 
per. However, we note that all the 18 systems in common 
between Planck and ROSAT have an average temperature 
higher than 5 keV, and thus in these systems the temper- 
ature is not expected to fall below ^ 2 keV. Conversely, 
this effect may affect the average profiles, which include a 
number of low-temperature systems (see E12). 



3.2. Temperature 

We can see in Fig. [3] that the combined 
temperature profiles extend out to ^ 2i?5oo, 
to the virial radius of these systems {Rvir 
in a ACDM cosmology), and a factor ^ 
yond the previo us results obtained f r om rel a tively large 
clust e r samples (iLeccardi fc Molendi 20081 j: IPra tt etal] 



X-ray/SZ 
i.e., out 
~ i?iio 

' 2 be- 



2001 IVikhlinin et al.M2006t iDe Grandi fc M olcndil l2002t 



Markevitch et all 1998[) . thus probing a volume ^ 8 times 
larger. The temperature is found to decline by a factor ~ 2.5 
from 0.2i?5oo to i?200 ^ l-52i?5oo, and to follow a power-law 
decrease like T{r) ^ r~", with a = 0.76 ±0.03 in the radial 
range [0.5 — 2]i?5oo. A slightly flatter proflle is obtained for 
the median of the 18 individual systems (a — 0.51 ± 0.05). 
This may indicate that different selection biases are affect- 
ing the Planck and ROSAT samples. In any case, the me- 
dian data points are found to agree rather well at all radii. 
Since the pressure and the density are obtained from the 
same systems, the value obtained by averaging the 18 indi- 
vidual objects is preferred. 

We also computed average temperature profiles for the 
CC/NCC populations separately, to look for differences be- 
tween clusters states. The resulting profiles are shown in 
Fig. m Both for the average profiles and the individual sys- 
tems, we observe no significant differences between the two 



Fig. 5. Average self-similar scaled entropy profiles obtained 
from the median of the 18 individual systems (black points) 
and the average pressure and density profiles (green shaded 
area). For comparison, the dashed line s hows the relatio n 
expected from gravitational collapse only (jVoit et al.ll2005[ ). 
using the proper normalization (Eq. [TU]). 



cluster populations. Beyond ^ 0.4i?5oo, NCC clusters are 
known to exhibit a higher density (by ~ 15 %) than CC clus- 
ters, which was highlighted in E12 (see also lMaughan et"aI1 
I2012i ). When the central slope 7 of the GNFW profile (Eq. 
[5]) is left free in the MCMC procedure (see Sect. 12. 21) . a simi- 
lar effect is seen for the pressure, such that the temperature 
profiles for the two populations are found to agree. 

3.3. Entropy 

In addition to the gas temperature, we used Eq. [5]to probe 
the behavior of the gas entropy at i?5oo and be yond. Several 
recen t results using the S u zaku s atellite (e.g.. iBautz et al.l 



2009"; 'Kawaharada et al.l 120101: ISimionescu et al.l 12011 
Walker et al. 2 012a ) found that the gas entropy increases 
less steeply than expected from gravitatio nal collapse {K 



„i.i 



, IVoit et all [20051 : iBorgani et aIll2005D. sugges ting that 



non-gravit ational effects (e.g ., iLapi et al.l I2O1O0 or gas 
clumping (jNagai fc Laul [20111 ) play an important role be- 
yond ^ i?5oo- Our dataset thus offers an excellent oppor- 
tunity to test the validity of these results in the frame- 
work of a large sample. To allow a meaningful comparison 
of our entropy profiles with the theoretical expectations, 
we rescaled them using the self-similar expectation (Eq. 
|4]) and compared them with the universal profile obtained 
from non-radiative simulations (jVoit et al.ll2005[ ). which is 
normalized by the quantity 



^^500 = 106 



M; 



500 



(i014Mg 



2/3 / -, N 2/3 

i—j h{z)-^/^ kcV cm^, (10) 



where h = flh/i^m. is the cosmic baryon fraction (see Eq. 3 
and 4 of lPratt et al.ll2010l for details). In Fig.[5]we show the 
average entropy profiles obtained from the universal pres- 
sure and density profiles, and from the median of the 18 
individual systems. The same figure showing the profiles ob- 
tained with the two deprojection techniques can be found in 
Fig. IC.5l of the online material. These pr ofiles are compare d 
with the self-similar expecta tion from Voit et al.l (l2005l). 
scaled according to Eq. [TU] ([Pratt et aL 2010t [Sun et al.[ 
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Fig. 6. Average entropy profiles for the CC (blue) and the 
NCC (red) populations. The shaded areas represent the 
quantities obtained from the average pressure and density 
profiles, while the data points show the median of the 6 
CC and 12 NCC individual systems. For comparison, the 
dashed line sh ows the relation expected from gravitational 
collapse only (|Voit et al.l[2005D . using the proper normal- 
ization (Eg.lTO)). 



1200911 . As can be seen in Fig. [5j the combined X-ray/SZ 
profiles are indeed fiatter than the self-similar relation in 
the radial range [0.5 — 2]i?5oo. Fitting the average entropy 
profile with a power law {K oc r"), we find a typical slope 
a = 0.60±0.03 in this radial range. As already found for the 
temperature, the individual systems are returning a some- 
what steeper slope, a = 0.90 ± 0.06. This may again reflect 
differences in the selection of the two samples. However, as 
can be see i i in F ig. [SJ the differences with respect to the 
IVoit et aTl (|200H) proflle tend to decrease when going far- 
ther out. This result disagr ees with the behavio r generally 
observed with Suzaku (see lWalker et al]l2012bl and refer- 
ences therein). 

In Fig. [S] we show the average entropy proflles for the 
CC and NCC populations, derived both from the average 
pressure and EM profiles and for the 18 individual objects. 
The profiles are compared to the self-similar expectation 
r^-^. The power-law slopes obtained in the [0.5 — 1] and 
[1 — 2]i?5oo ranges are presented in Table [TJ We notice a sig- 
nificant difference between the CC and NCC slopes in the 
outer regions. In CC clusters, the slope of the entropy pro- 
file is consistent with the self-similar expectations, whereas 
in NCC clusters we observe slightly flatter entropy proflles. 
Similar results are found for the average of the individual 
systems. This result is explained by the h igher density in 
the o uter regions of NCC systems (E12; iMaughan et al.l 
I2012D . 



3.4. Polytropic index 

Another quantity of interest for compa rison with semi- 
analytical models is the polytropic index (F "ozzi fc NormanI 
120011: lAscasibar et"all 120031: ICapelo et al.ll2012D . which i 
defined as 



IS 



dlogP 
dlogp ■ 



Table 1. Slope of the average entropy profiles in the ra- 
dial range [0.5 — 1] and [1 — 2]_R5oo obtained from different 
datasets. 



Dataset 



^[0.5--11 



Average 
18 objects 
Average CC 
6 objects CC 
Average NCC 



0.59 ±0.06 
0.86 ±0.11 
0.74 ± 0.06 
0.80 ±0.12 
0.58 ±0.08 



0.82 ±0.11 
1.01 ±0.21 
1.14 ± 0.16 
1.16 ± 0.24 
0.85 ±0.11 



12 objects NCC 0.74 ± 0.16 1.05 ± 0.30 



0. 
0. 



2 
1.5 
1 

0.5 


-0.5 
-1 
-1.5 




I 

-2.5 



log(n h{z) 



(11) 



Fig. 7. log P — log n diagram for the CC (blue triangles) 
and NCC (red circles) cluster populations. The solid lines 
show the best fit to the data using a rolling index. 



The polytropic relation that relates gas density and pres- 
sure by a simple power law is the simplest extension of the 
isothermal model and is indeed expected in ideal adiabatic 
processes. In these cases, the polytropic index is equal to 
5/3 and represents the ratio of specific heats at constant 
pressure and constant volume. As a result of the combined 
action of shock heating and adiabatic co mpression, the in- 
dex 7 is found in numerical studies fe.g.. lTozzi fc NormanI 
120011) to be between 0.8 and 1.2 over the cluster's volume. 
It has been shown (e.g., lAscasibar et al . 2003) that poly- 
tropic models provide a fairly accurate description of the 
radial structure of galaxy groups and clusters, failing in re- 
producing the cluster properties in merging systems, where 
the dark matter potential may differ considerably from a 
predicted NEW form and the assumptions of hydrostatic 
equilibrium might not hold. 

Since pressure and density are our primary observables, 
the dataset explored here is ideal to measure the polytropic 
index and study its dependence with radius. In Fig. |7| we 
show the median pressure as a function of the median gas 
density, for the 6 CC and 12 NCC systems independently. 
An orthogonal technique was used to fit the data taking the 
uncertainties in both axes into account. The best- fit values 
polytropic indices are given in Table |21 

We measure 7 = 1.21 ±0.03 (1.20 ±0.03) for CC (NCC) 
systems over the radial range 0.1 — 2i?5oo, with no significa- 
tive dependence upon the radial range considered. Overall, 
these values show a significant departure from the isother- 
mal condition in a systematic way at any radius, confirming 
an effective polytropic index of about 1.2 predicted as the 
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Table 2. Best- fit polytropic index for various datasets and 
radial ranges. 



Dataset 


70 


1- 


0.5 


70.5- 


-1 


71- 


2 


Entire sample 


1.06 


± 


0.07 


1.24 ± 


0.07 


1.24 ± 


0.11 


CC 


1.11 


± 


0.10 


1.22 ± 


0.14 


1.23 ± 


0.16 


NCC 


1.05 


± 


0.12 


1.26 ± 


0.12 


1.31 ± 


0.14 



outcome of the processes of shock heating and adiabatic 
compression acting on the X-ray emitting plasma during 
its accretion in the cluster's potential well. 



4. Discussion 

4.1. Comparison with previous results 



Exploiting the low background of the Suzaku satellite, sev- 
eral recent X-ray works claimed the mea surement of the 
ther modynamic properties out to i?200 (^see lAkamatsu et al.l 
1201 It IWalker et a"nl2012bl and references therein). These 
results are reporting a decline in the temperature out to 
R200, in qualitative agreement with the results presented 
here (see Sect. 13.21 and Fig. [3]), and a fiattening of the en- 
tropy profiles at the outermost radii, falling below the pre- 
diction from numeri cal simulations. For relaxed systems, 
IWalker et al.l (|2012b[ ) suggested the existence of a univer- 
sal entropy profile, showing a general flattening of the en- 
tropy profiles. Using a method similar to the one presented 
here, they combined the universal pressure profile from 
P12 with the mean deprojected density profile from E12, 
and claimed that these measurements are supporting the 
Suzaku measurements in the case of relaxed clusters. This 
claim is however incorrect. Indeed, IWalker et al.l (l2012b[ ) 
used the average profiles for the entire P12 and E12 sam- 
ples, and compared them with results obtained by Suzaku 
on relaxed systems. Our analysis (see Fig. [S] and M and 
Table [T]) clearly shows that the trend of entropy flatten- 
ing is driven by the NCC systems, which are in majority 
both in the Planck and ROSAT samples. The comparison 
with the relaxed Suzaku systems is therefore not appro- 
priate. Moreover, to compare with the simulation results, 
IWalker et a l. (2012b) arbitrarily set the normalization of 
the lVoit et al.1 ((2005) profile to match the observed profiles 
at a radius of 0.3i?2oo, while the normaliz ation (Eq.fTUj) can- 
not be fixed arbitrarily ()Sun et al.ll2"00l iPratt et al.lboiOt 
iHumphrev et al.ll2012D . This causes the CC profiles to fall 
below the self-similar prediction at larger radii, whereas our 
analysis shows that, in particular for CC clusters, the pro- 
files are in e xcess compared to the p redicted entr opy o ut to 
~ i?50o (see ICavagnolo et al.l [20091 iPratt et"alll2010l and 
references therein), but converge to the expected behavior 
in the outermost regions. 

4.2. Comparison witti ENZO simulations 

We compared our observational results with the proper- 
ties of simulated galaxy clusters produced in a set of high- 
resolution cosmological sim ulations with the ad aptive mesh 
refinement code ENZO (see IVazza et al.l l2010l for the de- 
scription of the simulated cluster sample). This set of 
simulations already allowed us to highlight some impor- 
tant differences in the amount of turbulent motions, shock 



waves, azimuthal scatter, and gas clumping factor between 
CC-like and NCC-like simulated systems, showing that all 
these effects tend to be s tronger in NCC-like clus te rs com- 
pared to CC-like o bjects (jVazza et al.ll2niflll2011alihl . l2ni2bl: 
lEckert et al.ll20l3 ). From the results of El 2, we report that 
the gas density profiles of CC-like and NCC-like clusters 
with simple non-radiative setup qualitatively agree with ob- 
servations for > i?5oo , while present results are at odds with 
observations in the internal radial range. Similar trends are 
also found in the simulated gas pressure. Overall, these re- 
sults may be understood in terms of a slightly different 
gas mass distribution, due to the lack of radiative cooling 
and inward contraction of the gas/DM atmosphere in non- 
radiative systems (e.g., E12 and discussion therein). This 
comparison can now be extended to all other thermody- 
namic properties of simulated and observed clusters out to 

For each cluster, we computed the profile of mass- 
weighted temperature, Tmw, and mass- weighted entropy, 
Tmw/n^^^, and generated average profiles for the CC-like 
and NCC-like classes at 2: = using the same weighting 
by the cluster mass, as in Eq. Figure [5] shows the av- 
erage temperature profiles obtained in this way compared 
with the observational results of the previous sections (me- 
dian of the 18 objects). Looking at the mass- weighted tem- 
perature for > 0.1i?5oo (where the effect of missing cool- 
ing in the simulations should be unimportant), the degree 
of qualitative agreement with observations is similar for 
both populations, with a tendency of simulated CC-like 
systems to present a slightly steeper trend than the ob- 
served CC. Apart from this, the agreement between the 
observed and simulated temperature profiles is excellent. 
In the right hand panel of Fig. [8] we show the average en- 
tropy profiles obtained in this way compared with the ob- 
servational results of the previous sections (median of the 
18 objects). While inside 0.2 i?200 the trend of simulated 
and observed NCC are qualitatively similar (with NCC sys- 
tems showing a flatter profile and a higher level of central 
entropy), the relative difference between CC and NCC is 
much more marked in real systems, which points toward 
the additional mechanisms of radiative cooling and non- 
gravitational heating. For i? > 0.2 i?2oo the slope of the 
simulated entropy profiles is close to the self-similar expec- 
tation, but with the tendency of producing a ^ 20 — 30 
percent higher entropy at all radii. This finding is consis- 
tent with the result we obtained in E12, where the gas den- 
sity at ~ i?200 in our simulated clusters has been shown to 
be slightly lower compared to the density profiles obtained 
with ROSAT. In addition, unlike in observations, the simu- 
lated NCC systems present no evidence of a flatter entropy 
profile at large radii. 

4.3. Entropy flattening or entropy excess? 

In the central regions, galaxy clusters are known to ex- 
hibit an excess of entropy with respect to the prediction 
from pure gravitational collapse, which indicates the exis- 
tence of an additional heating mechanism, such as feedback 
from active galactic nuclei (AGN) within already-for med 
clusters (see lFabianll2012t iMcNamara fc NulsenI l2012l for 
reviews ) or intense h eating episo des until the formation 
epoch (jMcCarthv et al . 2008; B urns et al.1 [2008). This ef- 
fect is known to be stronger in group -scale objects (e.g., 
iDavid et al.l Il996t iPonman et all I1999D , which causes the 
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Fig. 8. Comparison of average temperature (left) and entropy (right) profiles for the observed CC (cyan) and the NCC 
(yellow) populations derived from the median of the 18 objects, and the CC-like (blue) and NCC-like (red) populations 
simulated with ENZO. The dot-dashed lines show the median profiles from our simulated ENZO sample, while the solid 
lines indicate the confidence intervals on the median. In the right panel, we additionally show the self-similar expectation 
from lPratt et alJ (|20100 (dashed gray line). 



inner slope of the entropy profiles to be flatt er in low- 
mass systems with respect to massive clusters (jSun et al.l 
l2009t) . The entropy excess is known to extend si gnificantly 
beyon d the cluster core, at least out to i?25on (ISun et al.l 
200911 and possibly ev en out to i?5oo (jPratt et al. 2010t 



Humphrev et al.l[2012[ ). The results presented here for re- 



laxed systems indicate that instead of a flattening at large 
radii, the average entropy profiles actually match the self- 
similar expectation at ^ i?500: and are consistent with the 
expected behavior beyond this radius (see Table [1]) . The 
flatter slopes observed in the outer regions are thus result- 
ing from the entropy excess in the central regions and ex- 
tending out to large radii, instead of highlighting a further 
flattening at the outermost radii. 

This result is at odds with most of the results ob- 
tained from Suzaku data (e.g., Simionc scu et al.l 1201 it 
lAkamatsu et al.l 1201 ll : IWalker et al.l I2012al') . It highlights 
the importance of using the expected normalization for the 
self-similar entropy profile when comparing with the pre- 
dicted behavior in stead of normalizing at an arbitrary ra- 
dius, as shown by iHumphrev et al.l (j2012f) for the galaxy 
group RX J1159-t-5531, which exhibits a behavior very sim- 
ilar to our findings. Alternatively, the difference could also 
arise from the fact that most of the Suzaku measurements 
were obtained along narrow sectors, which may not be rep- 
resentative of the global properties in the virial region. 
Indeed, E12 demonstrated that deviations from azimuthal 
symmetry increase with radius, in particular in CC sys- 
tems. A coverage of the entire azimuth is thus required to 
infer the average properties of a cluster at large radii. In 
addition, we stress that, un like the case of XMM-Newton 
(|Leccardi fc Molendill2008l3 ). the level of systematic uncer- 
tainties in recovering the thermodynamical properties at 



large radii from Suzaku data is still poorly known, and no 
clear criteria are available to assess the reliability of a result 
(Molendi, private communication). All these effects could 
contribute to some extent to the discrepancy between the 
results presented here and those of Suzaku. However, deter- 
mining the exact origin of the discrepancy would require a 
detailed analysis of each case separately, which is beyond 
the scope of this paper. 

4.4. Correcting for the gas depletion 



In the REXCESS sample. IPratt et al.l (|2010[ ) observed that 
the excess entropy correlates with the gas depletion with 
respect to the expected value, such that when the entropy 
profiles are corrected by the observed gas fraction profile 
iifgas/fb)'^^^, where fb is the expected hot gas fraction), 
the expectation from gravitational collapse is recovered. In 
other terms, non-gravitational effects (AGN feedback, pre- 
heating) injects a substantial amount of entropy per parti- 
cle, which leads to a depletion of gas in the central regions 
and a high central entropy. Correcting for the gas deple- 
tion allows one to correct for the redistribution of the gas 
induced by non-gravitational effects. 

The dataset explored here offers the unique opportu- 
nity to test whether this effect also applies to the external 
regions. Using the gas fraction profiles presented in Paper 
II, we rescaled the rec onstructed entropy profile using the 
method introduced by iPratt et aD (|2010[ ) , both for the av- 
erage profile and the 18 individual systems. In Fig. we 
show the entropy profiles scaled by the self-similar 7^500 
(Eq. [To)) and correct e d for the variations in gas fraction 
following jPratt et al.l ()2010D . As we can see in the figure, 
the result obtained within i?5oo bv lPratt et al.l (|2010D holds 
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Fig. 9. Average self-similar scaled entropy profiles for the 
CC (blue) and NCC (red) cluster populations, corrected 
for the varia,tion o f fgas with respect to the expected value 
(jPratt et al.ll20ld[) . The dashed line shows the e ntropy pro- 
file e xpected from gravitational collapse only (jVoit et alj 
120051) . 



over the entire radial range (see also lHumphrev et al.ll20l3l , 
in particular in the case of relaxed CC clusters. For NCC 
clusters, slight deviations with respect to the prediction are 
observed, although the general trend is qualitatively sim- 
ilar. In these systems, we observe a gas fraction in excess 
of the expected hot gas value beyond i?5oo (see Paper II), 
which can be caused either by the breakdown of hydro- 
static equilibrium under the influence of non-gravitational 
effects (e.g., turbulence, bulk motions, cosmic rays), or 
by an inhomogeneous gas distribution (see Paper II and 
references therein). Correcting for this observed excess of 
baryons in cluster outskirts thus allows one to recover the 
expected behavior in the outer regions of perturbed sys- 
tems. Conversely, in CC systems both the observed hot gas 
fraction and the entropy (Fig ^ agree with the expected 
behavior, indicating that on the large scale these systems 
are governed mainly by simple gravitational physics. 

4.5. Possible interpretations 
4.5.1. Gas clumping 

As we can see in Fig. [Bl the entropy of NCC clusters falls 
below the expectation from pure gravitational collapse be- 
yond ~ 0.7i?5oo- Since the X-ray emissivity is proportional 
to the squared density, inhomogeneities in the gas distri- 
bution can lead to biased density measurements if one as- 
sumes a constant density within each spherical shell. This 
could bias the measured entropy low, which would lead to 
a d eficit of entropy with respect to the exp e cted behav- 
ior (iMathiesen et al.l[l999l: iNaeai k Laull201lllVazza et all 
l2012b t). The lower entropv of NCC systems could thus indi- 
cate the presence of significant gas clumping at large radii. 
In the same way, large deviations from azimuthal symme- 
try could affect the NCC population, although the mean 
azimuthal scatter on large scales appears to be similar in 
both cluster populati ons (see E12). 

Following Mathi esen et al.l (|1999D . we define the clump- 
ing factor as the bias in the density reconstructed from 
X-ray measurements compared to the true density, C = 



Assuming th at the diffe r ence b etween the NCC 
entropy profile and the IVoit et ahl (|2005h profile is caused 
by gas clumping, one can thus derive the clumping factor 
needed to recover the self-similar expectation. At i?200: we 
find VC = 1.23 ± 0.06 (1.18 ± 0.12 from the 12 individual 
systems). This value is in the range allowe d by simulations 
(liviathiesen et al.lll999 ': 'Nagai fc Lau"201lV and well below 
the value proposed by Simionescu et al.. (.201111 to explain 
the Suzaku observations of the NW arm of Perseus. This 
suggests that either this measurement, which covers less 
than 5% of the azimuth at i?200: is not representative of 
the cluster as a whole, or another mechanism must be in- 
voked to reconcile it with the theoretical expectations. 

Numerical simulations also predict differences in the 
clumping factor between merging systems and relaxed clus- 
ters, even if by a smaller factor than required here. Since 
merging clusters exhibit more accretion than relaxed sys- 
tems, especially along filaments, they arc expected to host 
a significantly large r number of clumps (Nagai & Lau 201l|; 
IVazza et al]|2012bl ). in agreement with the trend observed 
here for NCC. Conversely, in the case of CC systems good 
agreement is found with the predictions of gravitational col- 
lapse of uniform smooth cells of matter, which translates 
into a clumping factor consistent with unity at all radii. 



4.5.2. Non-thermal energy supply 

Cosmological simulations generally agree in reporting that 
the amount of turbulent energy is expected to be ~ 20 — 30 
percent of the thermal energy at i?200 iu NCC systems, and 
only ^ 15 percent in CC (e.g.. lLau et al.ll2009t IVazza et al.l 
l2011al ). After major merging events, the magnetic energy 
in the ICM should b e boosted by a factor ^2 — 4 (e.g ., 
iRoettiger etHI 119991: iDolag et all l2002t IXu et al.l 120111 ). 
and so is the relative boost of CR energy within clusters 
following merger sh ock waves (e.g., iPfrommer et al.l 120071 : 
IVazza et al.ll2012al) . Indeed, the more intense activity of 
merger and large-scale accretions (e.g. from filaments) that 
is likely at work in NCC systems is expected to promote a 
larger injection of turbulent motions, cosmic-ray (CR) en- 
ergy and magnetic field amplification compared to CC sys- 
tems, soon after crossing the cluster virial radius. All these 
mechanisms reduce the thermalization efficiency of kinetic 
energy within clusters so that, for a given total mass, an 
NCC system is expected to have a lower temperature than 
a CC system in the outskirts (< R200), where non-thermal 
phenomena are expected to be more important due to the 
departure from the simple "smooth accretion" at work in 
real s ystems. This is indeed what variou s simulations re- 
port ( Lau et aT1l2009t fVazza et al.ll2011a^ . including those 
presented in this paper (see Sect. 14. 2p . 

Conversely, our observations show that, outside the 
core, the CC and NCC temperature profiles are indistin- 
guishable (see Fig. 2]). This would seem to imply that the 
thermalization process in the outer regions of NCC and 
CC systems is equally efficient, possibly indicating that the 
mechanism that are responsible for thermalization are sig- 
nificantly more rapid than those currently included in sim- 
ulations. While this is an intriguing possibility, we must 
however note that, given the limited statistical quality of 
our data, we cannot exclude that NCC profiles are on aver- 
age lower than CC profiles by about 20-25% at i?200- Thus, 
while our data do not require that NCC and CC system 
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have a different fraction of non-thermal energy they can- 
not, at least to some extent, exclude this. 

5. Conclusion 

Combining X-ray (ROSAT) and SZ (Planck) data, we an- 
alyzed the properties of the ICM at large radii using both 
sample-averaged profiles and the data for 18 individual sys- 
tems. This allowed us to measure for the first time the ther- 
modynamic properties (temperature, entropy) of galaxy 
clusters out to the virial radius in relatively large samples. 
Gravitating mass profiles obtained using the same dataset 
are presented in a companion paper. Our results can be 
summarized as follows: 

— We have shown that the thermodynamic properties of 
the ICM can be efficiently recovered by using the ther- 
mal pressure profiles from SZ observations and the gas 
density profiles from X-ray data. Comparing the results 
with temperature profiles from the literature, we found 
that this method allows us to reconstruct the gas tem- 
perature within ~ 10% of the X-ray spectroscopic value. 

— In cluster outskirts, the temperature is found to de- 
crease like T(r) c>c 7'-o.5i±o.06 |-j^g radial range 
[0.5 - 2]i?5oo (r(r) oc r~°-''^^°-°^ from the average pro- 
files). The temperature at i?200 is on average lower than 
that at 0.2i?5oo (which is close to the spectroscopic X- 
ray temperature) by a factor of 2.5. No significant dif- 
ferences are observed between the temperature profiles 
of relaxed CC clusters and unrelaxed NCC objects. 

— As indicated by several recent Suzaku observations, the 

— 5 /3 

entropy K = Prie is slightly flatter than the self- 
similar expectation in the radial range [0.5 — 2]i?5oo, 
although it rises steadily with radius. However, we 
note that at the outermost radii the observed behav- 
ior of the entropy is close to the one expected from 
simulations, once th e proper normalization for the en - 
tropy profile is used (I Voit et al.ll200l iPratt et al.ll2010tl . 
When a correctio n for the gas depletion is introduced 
(jPratt et al.ll2010f ). the entropy profiles for the two pop- 
ulations are found to be in close agreement with the ex- 
pectation from gravitational collapse. Therefore, instead 
of a decrement of ent ropy at la rge radii as indicated by 
Suzaku (see also Walker et al.l [2012b'). our results indi- 
cate that the well-known entropy excess observed in the 
inner regions extends well beyond the core of clusters 

out to i?200- 

— Unlike what was recently reported by I Walker et al.l 
(2012b), the average entropy profile in relaxed systems 
is found to agree with the expectations from gravita- 
tional collapse at the outermost radii, while in per- 
turbed NCC systems the average entropy profile falls be- 
low the self-similar prediction. This indicates that non- 
gravitational effects (turbulence, cosmic rays) and/or 
accretion of non-virialized structures (i.e., gas clump- 
ing) is more important in these systems. 

— The polytropic index is found to exhibit a weak radial 
dependence and to have an average value of 1.2, with no 
obvious difference between the two cluster populations. 
This result agrees well with the theoretical expec tations 
fe.g.. l^zzi fc No rma ^ I20nit ICanelo et al.ll2nil . 

— The differences observed between CC and NCC clusters 
at large radii could be explained by a higher clump- 
ing factor in the latter population. If the deviations 



with respect to the self-similar expectations are due 
to gas clumping, we measure a clumping factor VC = 
1.23±0.06 at i?200 for NCC systems, which is well within 
the range expected from numerical simulations. In CC 
systems, our observations are consistent with a smooth 
ICM at all radii. 
— Alternatively, some of the observed trends in CC and 
NCC profiles could be interpreted in terms of a more ef- 
ficient injection of non-thermal energy (turbulence, cos- 
mic rays, magnetic fields) in the latter population fol- 
lowing from merging events, at the level of ^ 20 — 40 
percent at i?2oo- However, given the lack of observed dif- 
ferences in temperature between CC and NCC systems, 
this interpretation is not favored by our data, although 
it cannot be firmly excluded. 
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Appendix A: Master table 

In Table ETTl we give the list of clusters in common between 
the Planck and ROSAT sampl es. For the detailed proper- 
ties of these clusters, we refer to lPlanck Collaboration et al.l 
(|201lh and E12. 




2^:1 



Fig. B.l. Planck / ROSAT reconstructed temperature as a 
function of the obser ved temperat ure from XMM -Newton 
(jSnowden et al.ll2008[ ) or Chandra (jCavagnolo et al., ,20091) 
for the 18 objects in common between the two samples. The 
dashed red curve represents the best-fit function Tree = 
0.97Tobs. 

Appendix B: Comparison with literature 
temperature profiles 

To ensure that the reconstruction of thermodynamic quan- 
tities from the X-ray density and the SZ pressure is giving 
correct results, we used the data for the 18 individual sys- 
tems in common between the ROSAT and Planck samples 
(see Sect. 12.1 1 and Table ETl]) and com pared them with tem- 
perature pro files from XMM-Newto n ( Snowden et al.l[2008l) 
or Chandra (jCavagnolo et al.ll2009D . Given that the angu- 
lar resolution of ROSAT is poorer than that of Chandra 
and XMM-Newton, we restricted our comparison to the ra- 
dial range beyond 0.1i?5oo- We also excluded the observa- 
tional data points beyond i?5oo, since the techniques used 
to extract X-ray temperatures in large samples may not be 
valid in regions with a low signal-to-background ratio (see 
Ettori & Molendi 2011). 

In Fig. lB.ll wc show the temperatures reconstructed us- 
ing our method versus the observed temperatures for all 
available data points. The data were fit with a maximum- 
likelihood estimator that accounts for the scatter around 
the mean value, which needs to be computed simultane- 
ously (see lMaccacaro et al.ll988l and paper II). The best-fit 
value for the data points in the [0.1 — l]i?5oo range gives 



T 



To, 



= 0.97 ±0.01, 



6s 



(B.l) 



with a scatter of 9%. This analysis demonstrates that com- 
bining X-ray density profiles with SZ pressure profiles is an 
effective way of reconstructing thermodynamic quantities. 

Appendix C: Notes on individual objects 

— A 478: The line of sight of this system is known 
to exhibit a complex, varying column density 
(|Pointecouteau et"aLl l2Q04f ). For this work, as in 
E12, we fixed Nh to the 21cm value (jKalberla et al.l 
|2005(). This may result in a slightly incorrect density 
profile. 
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Table A.l. Individual systems in common between the E12 and P12 samples. 



Name 


State 


qt 




OLK 




fgas,500 




/gas, 200 


R, 


max 


A85 


CC 


0.05 ±0.08 


1, 


.52 ±0. 


,09 


11.3 ± 1, 


,1 


9.4 ± 1.0 


1, 


,59 


A119 


NCC 


-0.41 ± 0.13 


0, 


.94 ±0. 


,14 


15.3 ±2, 


,4 


13.9 ± 2.5 


1, 


,55 


A401 


NCC 


"0.27 ±0.08 


1, 


.02 ±0. 


,08 


15.8 ± 1, 


,7 


18.9 ± 2.1 


1, 


,83 


A478 


CC 


-0.40 ±0.13 


0, 


,94 ± 0, 


13 


13.4 ± 1, 


,7 


15.7 ±2.3 


1, 


,81 


A665 


NCC 


-0.54 ±0.10 


0, 


,86 ± 0, 


,11 


20.7 ± 2, 


,6 


24.8 ±3.4 


2, 


,11 


A1651 


NCC 


-0.09 ±0.10 


1, 


,35 ±0. 


,10 


12.8 ± 1, 


,3 


11.8 ± 1.5 


1, 


,43 


A1689 


NCC 


-0.66 ±0.11 


0, 


,72 ± 0. 


,12 


13.9 ± 1, 


,9 


20.9 ±3.3 


2, 


,10 


A 1795 


CC 


-0.60 ±0.14 


0, 


,95 ± 0. 


,15 


10.9 ± 1, 


,4 


12.3 ± 2.2 


1, 


,56 


A2029 


CC 


-0.44 ± 0.09 


0, 


,97 ± 0. 


,10 


13.8 ± 1, 


,5 


15.8 ±1.8 


1, 


,62 


A2163 


NCC 


-0.29 ±0.09 


1, 


,15 ±0. 


,11 


32.5 ± 3, 


,0 


27.3 ±4.0 


1, 


,53 


A2204 


CC 


-0.61 ± 0.09 


0, 


,70 ± 0. 


11 


14.9 ± 1, 


,7 


20.1 ± 3.2 


1, 


,80 


A2218 


NCC 


-0.16 ±0.10 


1, 


,10 ±0. 


,10 


19.6 ± 1, 


,6 


22.6 ±2.5 


2 


,07 


A2255 


NCC 


-0.54 ±0.10 


0, 


,72 ± 0. 


,10 


11.7± 1, 


,4 


15.5 ±2.3 


2, 


,11 


A2256 


NCC 


-1.02 ±0.13 


0, 


,42 ± 0. 


,10 


14.2 ± 3, 


,0 


21.7 ±4.4 


1, 


,81 


A3112 


CC 


-0.80 ±0.10 


0, 


,45 ± 0. 


,10 


11.8 ± 1, 


,3 


16.5 ±2.3 


1, 


,75 


A3158 


NCC 


-0.57 ±0.12 


0, 


,77 ± 0. 


,13 


17.5 ±2, 


,6 


22.8 ±4.0 


1, 


,68 


A3266 


NCC 


-0.22 ±0.12 


1, 


,64 ± 0. 


,14 


12.3 ± 1, 


,7 


10.0 ±1.9 


1, 


,51 


A3558 


NCC 


-0.14 ±0.10 


1, 


,34 ± 0. 


,11 


13.9 ±1, 


,7 


12.3 ± 1.6 


1, 


,66 



Column description. 1. Core state: cool-core (CC) or non-cool core (NCC). 2. Slope of the temperature profile in the radial range 
[0.5 — 1.5]-R5oo. 3. Same for the entropy. 4. Gas fraction at -R500 (see Paper II). 5. Same at ii200 . 6 . 2a maximum detection radius 
in the ROSAT data. 

— A4.OI: This system is in an early stage of merging with 
its neighbor A399, located ~ 35' SW, and connected 
to it by a bridge of X-ray emission. The density profile 
was thus extracted in a sector excluding the direction 
connecting A401 to A399. Because it is unclear whether 
the same was done for the SZ profile, the resulting 
profiles may be unreliable. 

— A2163: This very perturbed merging cluster may 
be out of hydrostatic equilibrium, which causes the 
mass estimates from X-rays to be overestimated by 
a factor of ~ 2 with respect to the mass derived 
from weak lensing and galaxy distribution (jSoucaill 
|2012[ ). Therefore, the scaling factors (Eq. H] and [T]) 
may be in correct, as well a. s the e stimate of i?5oo from 
iPlanck Collaboration et al.l (|20lH) . 
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Fig. C.l. Best fit to the Planck pressure profiles of the 18 systems (black points, from P12) with a GNFW profile (green 
lines), in the radial range of interest for this paper ([0.2 — 1.9]i?5oo)- The bottom panels show the ratio between the data 
and the model. 
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Fig. C.2. Same as Fig. [Cl 
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Fig. C.3. ROSAT density profiles of the 18 systems obtained through geometrical deprojection (black points) and 
parametric forward fitting (green-shaded areas), in the radial range of interest for this paper ([0.2 — 1.9]i?5oo). The 
bottom panels show the ratio between the two methods. 
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Fig. C.4. Same as Fig. [031 
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Fig. C.5. Average entropy profiles for the total population for the average populations and for the median of the 18 
systems in common. Similar to Fig. [5l but showing the differences between the parametric and deprojection techniques 
in both cases. 



